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We study the evolution of the fractal dimensions during anisotropic aggregation and disaggregation 
(i.e., when the applied field is switched off) processes, in a magneto-rheological fluid composed by 
super-paramagnetic particles with a diameter of f micron, suspended in water under the action of a 
constant uniaxial magnetic field. We use video-microscopy and image analysis to retrieve the form of 
the aggregates and calculate the following fractal dimensions: one-dimensional D\ , two-dimensional 
D2, perimeter-based fractal dimension D p and two-dimensional capacity dimension or box-counting 
dimension Db(2D). We apply this methodology for different values of the magnetic field amplitude 
and particle concentration. During aggregation, we can calculate an average value (D p } ~ 1.84 for 
the rod-like chains and we find a linear dependence of the capacity dimension with the ratio of two 
characteristic length scales. By means of a quantification of the roughness, we interpret D p as a 
measure of the clusters border roughness, while the capacity dimension represents how the chains 
fill the space. We also show how the quadratic deviation of the cluster contour height follows power- 
law behaviours in the initial stages of aggregation and disaggregation. We compare our results on 
projected fractal dimensions with previous theoretical and experimental studies about aggregation 
of magnetic particles. As a result, we check some recently proposed theoretical relation between D p 
and the three-dimensional fractal dimension Df. 

PACS numbers: 61.43.Hv, 83.80.Gv, 82.70.-y 



I. INTRODUCTION. 

The study of fractal geometry [l|, Q aiming to char- 
acterize the morphology of aggregates has been a wide 
and diverse research field on mathematics, chemistry, 
physics and biosciences along the last decades. For ex- 
ample, the concepts of fractal geometry have been useful 
in fields as diverse as the study of cancer growth and 
cell structure 0, 13] an d the distribution of galaxies @ . 
In the last three decades, several experimental studies 
about colloidal aggregation have been carried out aim- 
ing to determine fractal dimensions of colloidal aggre- 
gates, in special, by means of light scattering techniques 
[1) 0) The interest to determine fractal dimension 
of aggregates is related to the practical applications. In 
fact, by means of fractal dimension analysis, dynamic 
or morphological properties of aggregates can be investi- 
gated in processes relevant for water depuration, biology, 
medicine or engineering. For instance, properties such as 
deposition, sedimentation , settling [ll], [TJJ , rough- 



ness [l3( | . coagulation or porosity (see references on [l4l|). 
have been analysed in terms of fractal concepts. Two 
main classes of optical techniques are widely used in the 
study of the morphology of the ag greg ates in colloid sci- 
ence: light scattering techniques [151 ] and direct imag- 
ing techniques (e.g., video- microscopy) [ll]. Hereafter, 
we will study 2D-projected fractal dimensions by means 
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of image analysis in an experimental colloidal system, 
namely, a magneto-rheological fluid (MRF). We use an 
aqueous suspension of micron-sized super-paramagnetic 
particles that, under the action of a constant uniaxial 
magnetic field, acquire a magnetic dipole moment that 
forces them to aggregate into linear chains [T^llll Il9l. |2C| . 
Magnetic particles [2l|, [22|, [23| , magneto-rheological flu- 
ids [HI, [H, HI] , organic-solvent-based magnetic fluids [27| 
or magnetic liposomes [28| have been the object of sev- 
eral experimental investigations on fractal aggregation 
[29L I30I ] and fractal geometry. Studies on fractal dimen- 
sion by means of simulations have been also developed on 
magnetic particles (U [3l| or two-dimensional magnetic 
fluids simulations including van der Waals and electro- 
static forces [32j. In this paper, we use 2D images of the 
colloidal aggregates in a MRF to compute 2D-projected 
values of the fractal dimension of the aggregates. These 
projected fractal dimension gives a lot of information 
concerning the morphology and spatial properties of the 
clusters. Regarding fractal properties, many of the ex- 
perimental works for magnetic particle aggregation com- 
mented before, give different information: light scatter- 
ing techniques yield fractal dimension values related to 
the three-dimensional (3D) structure of the aggregate dis- 
tribution, while direct imaging techniques usually deal 
with two-dimensional (2D) projected images of the real 
three-dimensional structure of the aggregate distribution. 
However, some progress in this direction has been re- 
cently done [HI , and several possible relations have been 
theoretically found [HI, [H, [36| . We perform a detailed 
analysis of these studies on Section IVII 

The structure of the paper is as follows: after the In- 



2 




troduction, Sectionllll contains a description of the exper- 
imental setup and procedure. In Section HTT1 we include a 
brief recall on the definitions of the several fractal dimen- 
sions that will be used throughout the paper, including 
practical details about the computational methods. Sec- 
tion IIVI contains the description of the results and Sec- 
tion [V] displays some conclusions and prospects for future 
work. 



II. MATERIALS AND METHODS 

A. Experimental Setup. 

Our magneto-rheological sample consists on a com- 
mercial aqueous suspension of super-paramagnetic par- 
ticles, (Estapor Ml-070/60). These particles are com- 
posed by a polystyrene matrix with embedded magnetite 
crystals {Fe^O^) of small diameter (10 nm). Since these 
iron oxide grains are randomly oriented inside the micro- 
particles, the resulting average magnetic moment is neg- 
ligible in the absence of an external magnetic field. The 
super-paramagnetic particles have a radius a — 0.48 /im 
and a density p ~ 1.85 g/cm 3 with a magnetic content 
of 54.7%. 

The particles are denser than water and, therefore, 
they sediment to the bottom of the container. The height 
probability distribution of the particles can be calculated 
[33] and, for the particles here used, it can be shown that 
80 % of the particles lie in a narrow layer above the bot- 
tom glass plate, with a thickness of, approximately, twice 
the particle diameter. Hence, most of the particles are 
captured with a narrow field depth optical system, and 
the images correspond to an effective two-dimensional 
system. 

The surface of the latex micro-spheres is functional- 
ized with carboxylic groups. The suspension has been 



also provided with a 1 g/1 concentration of a surfactant 
(sodium dodecil sulfate). For these kind of suspensions 
we have measured a zeta potential by means of a Malvern 
Zetasizer Nano, obtaining values in between —110 mV 
and —60 mV for suspension pH in the range 6 — 7, re- 
flecting the presence of an electrostatic repulsive force in 
the system. 

The suspension sample is confined in a cylindrical cell, 
with a diameter of 6.5 mm and a height of 100 /Ltm, made 
of two quartz windows separated by a teflon spacer. A 
water jacket connected to a circulation thermostatic bath 
keeps the temperature constant during the experiments 
to T = 282 ± 1 K. The system used to generate the 
magnetic field is composed by a pair of coils fed by a 
power amplifier (BOP-50-4M) driven by a programmable 
arbitrary function generator (HP-33120A). This system 
allows for the generation of an uniform magnetic field 
along the measurement region (see more details about 
the experimental system on [38]). The video- microscopy 
setup used to record the images is located below the 
sample, and is formed by a long working distance mi- 
croscope (Navitar) with zoom capabilities attached to a 
CCD camera (Retiga-Ex) with 12 bit intensity resolu- 
tion and 1360 x 1036 px spatial resolution. The required 
illumination is generated from the top using an optical 
fiber Fostec light source. Full frame images are captured 
every 0.4 seconds during a total time of approximately 
8000 seconds at maximum spatial resolution. 



B. Experimental procedure. 

The control parameters in the experiments are the 
volume fraction of the suspension and suitably defined 
dimensionless magnetic interaction energy. Actually, 
when an external magnetic field H is applied, the par- 
ticles acquire a magnetic dipole moment aligned with 
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the external field direction. This magnetic moment is 
m = (4-7r/3)a 3 M, where M — xH, being M the mag- 
netization of the particle, and \ the particle's magnetic 
susceptibility. The relevant parameter in the study of 
these suspensions is the ratio between the magnetic in- 
teraction energy for two particles in contact, W m and the 
energy of thermal fluctuations, k B T: 

, _ W m _ /x m 2 
~ k B T ~ 16na 3 k B T [ ' 

where po is the vacuum magnetic permeability, k B the 
Boltzmann constant, and T the temperature. The val- 
ues of A in the experiment have been calculated from 
the magnetization curve of the particles that was ob- 
tained by means of a vibrating sample magnetometer. 
The value obtained for the saturation magnetization was 
46.5 emu/g. 

The volume fraction of particles <p m the solution is 
defined as the fraction of the volume occupied for the 
particles over the total volume of the solution. For the 
purpose of this study, it is more useful to use an effective 
surface fraction which is calculated by dividing the sum 
of the area of all cluster contained in the image by the 
image total area. Hereafter, (j>2D will refer to this effective 
surface fraction. 

These parameters, A and <p2D, can be used to define 
two characteristic length scales. First, we can define a 
distance, R\, at which the dipole-dipole interaction en- 
ergy is equal to the energy of thermal fluctuations, i.e., 
i?i = 2CA 1 / 3 . Second, we may define an initial average 

1/2 

inter-particle distance, Rq = 2a(\>^ D . The ratio of these 
two length scales allows us to distinguish between diffu- 
sion limited and field driven aggregation processes. If at 
the time the field is switched on R\ < Ro, the aggrega- 
tion process should be diffusion limited, while if R\ > Ro, 
the aggregation process should be field driven. 

The sedimentation velocity for the particles in the sus- 
pension can be calculated from the buoyancy force as 
v « (2a 2 g(p ~ p s ))/(9rj), where g is the gravity constant, 
and r\ and p s are the viscosity and density of the solvent, 
respectively. Thus the sedimentation velocity in our ex- 
periment is 0.43 fim /s. 

When the sample cell containing the magnetic suspen- 
sion is introduced in the experimental setup, it is left 
at least 15 minutes to equilibrate without external field, 
to allow for particle sedimentation and temperature sta- 
bilization. After this time, most of the particles lay in 
a thin layer of 2 particles diameter depth, located right 
above the bottom quartz window. 

After this equilibration time, we capture images dur- 
ing 5 minutes without field and then, a constant uniform 
magnetic field is switched on. The field triggers the ag- 
gregation process and, images are recorded during the 
time the field is being applied (approximately, 5000 sec- 
onds). Then, we switch off the magnetic field and capture 
another hour of images during the disaggregation pro- 
cess of the previously formed chains. Hence, we record 
approximately 21000 images for each experiment. 



1,0-1 l-H 1 • 1 1 1 1 1 1 1 . 1 1 1 , 16 



0,5- 



(a) i (b) ^ 

J J**" 


m •_ 

f -=C>(t) 











1000 2000 3000 4000 5000 6000 7000 0000 
time (sec) 

FIG. 2: Temporal evolution of mean cluster circularity (C) (t) 
and mean cluster size (S) (t) for a case of A = 1688 and 
<f>2D = 0.088. The magnetic field is applied at t — 300 s 
((a) region without field and (b) region with field where oc- 
curs the aggregation process) and disconnected at t = 4620 s 
((c) region without field and disaggregation process). 

This procedure is performed in experiments carried out 
at different values of A and 4>2D- The image analysis, data 
extraction and statistical calculations have been carried 
out with image processing software developed at our lab, 
based in the Image J[3!| free open source Java program 
for image analysis. In this software, we have implemented 
an adaptative threshold image processing algorithm that 
removes the image background and suitably captures the 
contour of the clusters. An example of a filtered image 
can be seen in FigfT] The contour of each cluster is then 
used to calculate geometrical properties such as cluster 
perimeter, P, area, A, Feret's diameter (longest distance 
between cluster contour points), If, and circularity, C — 
4ttA/P 2 . 

Moreover, we have calculated the cluster size proba- 
bility density, n s (t), which is the number of clusters of 
size s, per unit area, present in the system at time t [40l |. 
Then, n s (t), is used to calculate the average cluster size 
S(t) = J2s 2 n s (t)/J2sn s (t). 



III. FRACTAL DIMENSION USING IMAGE 
ANALYSIS. 

A fractal object presents self-similarity, i.e., when we 
change the scale where the object is observed, the result 
is very similar or identical to the original observed one. 
It is possible to classify fractals according to the varia- 
tions of this property. From a rigorous point of view, a 
fractal can be defined as a set of points whose Hauss- 
dorff dimension is larger that its topological dimension. 
In general, if we have a fractal quantity F, it is connected 
with a length / using a relation such as F ~ l Dc , where 
Dq is the fractal dimension associated with this fractal 
quantity. By means of fractal dimension, D f (also called 
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FIG. 3: Log-log plots of number of counts versus the in- 
verse of the box size for a box-counting method. Box sizes 
on the log-log plots are powers of 2, so the points on the 
graphs are equally spaced. We use a chain during aggregation 
(t = 4989.6 s) and during disaggregation (10 s after the field 
has been switched off) for an experiment with A =703 and 
(j> 2D =0.106. In the first case we obtain D B {2D) = 1.11±0.02 
and D B (2D) = 1.28 ± 0.02 in the second. 



This method gives the box-counting dimension or ca- 
pacity dimension (Dg hereafter). Classical box-counting 
method presents limitations that have to be taken into 
account: this method is very sensitive to the resolution 
and orientation of the image and the power-law relation- 
ship must be verified in a reasonable range of length 
scales to be able to state that fractal structure appears 
PH - An example can be seen on Fig. [3J where we cal- 
culate Db(2D) for a chain during aggregation and 10 s 
after the field has been turned off for an experiment with 
A =703 and 4>2D =0.106. In the first case, a capacity 
dimension Db{2D) = 1.11 ± 0.02 is obtained, whereas 
in the second, we obtain D B (2D) = 1.28 ± 0.02. The 
range where the linear regression is fitted is almost two 
orders of magnitude. In Section [IV! we use this method 
for several chains for each time and we calculate average 
values. 

Another method for obtaining fractal dimensions from 
these 2D images of the projected aggregates consists on 
the analysis of several projected fractal dimensions using 
the scale relation for fractal quantities explained before. 
We can consider various types of these fractal dimensions: 
one and two dimensional fractal dimensions, D\ and D2, 
which are obtained through the following form: 



.D3 or D), an interpretation can be made on how the ag- 
gregate fills the space or how primary particles fill the 
space of the volume occupied by the aggregate. An usual 
method to determine the fractal dimension of an aggre- 
gate knowing its number of particles, N , is through the 
radius of gyration, R g [29(. Then, the following expres- 
sion is verified: 



R g {N) ~ N 1/D f 



(2) 



This is a very reliable method for determining Df, if 
the number of particles for cluster is known with enough 
precision, being widely used in colloidal aggregation [4l|. 
However, it is usual in image analysis not to have a re- 
liable experimental setup or methodology for the right 
extraction of the number the particles on a cluster. In 
our case, we have searched for a compromise between 
reasonably good statistics and spatial resolution. As a 
result, we do not have enough image definition for de- 
tecting the individual particles inside each cluster (see 
FigOJ and only the contours of the clusters have been 
detected. 

An estimation method of Haussdorff dimension is the 
box-counting method consisting in covering the fractal 
object with boxes of different size and extracting the frac- 
tal dimension through a power law between the number 
of boxes and its size (42J. That is, if we cover a fractal 
objet with a number of boxes N(r) with side r, we can 
obtain its capacity dimension using the following expres- 
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lim 

r— >0 



log N(r) 
logr 



(3) 



I Dl - 



(4) 



where P represents the perimeter of a cluster, A the area 
and If the longest distance between two points in the 



cluster, the so-called Feret's diameter [34j. Another use- 
ful expression uses the area and the perimeter for obtain- 
ing the so-called perimeter-based fractal dimension D p in 
the following form: 



A ~ P 2 /°p 



(5) 



These relations let us calculate the fractal dimensions 
easily using 2D images captured with a digital camera. 
Thus, when the contour of the clusters or aggregates is 
detected by means of image analysis, we only need to 
calculate the perimeter, the area and the Feret's diame- 
ter for obtaining the corresponding 2D fractal dimension. 
However, some critics to the applicability of the area- 
perimeter method on Ec.([5]) have been lately reported 
[44J |. The main argument of that critic lies on the fact 
that digitalising can change the perimeter and area of the 
objects of study, if it is applied for objects oriented in 
different directions. In our case, this argument is not ap- 
plicable, because of the anisotropy of the clusters caused 
by the external magnetic field. An example of the calcu- 
lation of Dp using the area-perimeter method can be seen 
on Fig[4] The power law behaviour is verified in almost 
four orders of magnitude, showing that this method for 
obtaining fractal properties is adequate for a system like 
ours. 
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it is possible to connect the roughness exponent with a 
fractal dimension d |46i: 
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FIG. 4: Linear regression for the calculation of D p corre- 
sponding to a magneto-rheological fluid with concentration 
<f>2D = 0.088 and A = f688 after applying the magnetic field 
during 4500 s. A value of D p = 1.84 is obtained. 



A. Roughness. 

For testing the possible interpretation of the projected 
fractal dimensions as a reflection of the roughness of 
the boundary of the clusters, we compute the height 
of the clusters contour measured from a central line 
which crosses the clusters end to end which we named as 
hj(x,t). For each cluster j in a single image, this quan- 
tity depends of the position x of the contour point and 
time t. For this study, we calculate the average value, 
H{t) = (h(t))j, of this height for all the clusters con- 
tained at the corresponding image. This border rough- 
ness study is different to others previously made on a 
single magnetic chain [45| , where the fluctuations on the 
form of the chain were analysed, not its border. 

We also calculate the border height fluctuation by 
means of the interface width for the cluster j as: 



\ 



1 N 



(h(t))] 



(6) 



where i refers to the corresponding contour point and 
N is the total number of contour points. This quantity 
characterizes the roughness of the contour. We calculate 
a cluster average value of the border height fluctuation 
for each time as W(t) = (w(t))j. Also, we can define the 
height-height correlation: 



c(l) = d[h(x)-{h(x>))Y 



(J) 



where I = \x — x'\ is the distance between two contour 
points projections, x and x' , made into a central end to 
end line that crosses the cluster. This magnitude scales 
in the form c(l) ~ l a , being a the roughness exponent. 
This property has been widely used, for instance, for the 
study of surface growth (4(| and interfacial growth phe- 
nomena by electrochemical deposition [47j | . Moreover, 



d 



a 



(8) 



Later in this work, we will show as D p , the perimeter- 
based dimension, may be identified with this fractal di- 
mension d, confirming the morphologic nature of this 
quantity. We also use the temporal evolution of the av- 
erage height-height correlation Wit) for observing the 
changes on the cluster roughness during the aggregation 
and disaggregation processes. 



IV. RESULTS. 

We calculate the capacity dimension Db(2D) using 
several clusters in some images during aggregation and 
disaggregation. For an adequate calculation of this di- 
mension, we extract manually well-defined and represen- 
tative clusters in each image and calculate individually 
the box-counting dimension. Next, we make an average 
value of Db{2D) for each analysed time. This method 
is only valid when we study objects with enough size 
for fitting a linear regression in an acceptable range. Be- 
sides, because of the image resolution, we cannot use this 
methodology in the stages in which the number of free 
particles is predominant. We also obtain one-dimensional 
(Di), two-dimensional (-D2) and perimeter-based (D p ) 
fractal dimensions using Ec.(|4]) and Ec.©. Therefore, 
we fit three linear regressions in every captured image 
(every 0.4 sec) using our software, aiming to observe the 
changes on these fractal dimensions during aggregation 
and disaggregation. We require a regression correlation 
factor of at least r > 0.98 for considering that the result 
of its respective image is correct. This condition widely 
reduces the number of images that are useful for this 
study, specially when chains are not yet formed, that is 
when the particles are moving free in the liquid. As it 
can be seen on the example in Fig. [IJ we use a wide vi- 
sion field for detecting as many chains as possible. This 
is just what happens with Ds(2D), for small objects, 
like free particles, we do not have enough precision -or 
number of pixels- to observe clearly the particles' border, 
so we cannot suitably obtain their corresponding fractal 
dimension. 



A. Aggregation. 

As it has been commented in the Introduction, a 
few studies about fractal dimension of aggregates us- 
ing magnetic particles have been developed in the last 
two decades. The differences between experimental sys- 
tems imply that a complicated issue can arise from the 
comparison between different obtained fractal dimen- 
sions. However, it is obtained that all the values relat- 
ing chain-like structures of magnetic particles for three- 
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FIG. 5: Temporal evolution of projected projected fractal di- 
mensions during aggregation. Continuous lines D2 (bottom) 
and D p (top) for a case of A = 703 and <f>2D = 0.106 (cor- 
relation factor r > 0.98). The single points represent box- 
counting dimension Db(2D): A =131 and <j>2D =0.071 for 
squares; A =176 and 02D =0.115 for circles; A =703 and 
4>2D =0.106 for triangles. 



dimensional fractal dimensions are contained in a in- 
terval Df = 1.1 — 1.3. For instance, using light scat- 
tering, it is possible to obtain Df adjusting the scat- 
tering intensity / with the scattering vector q between 
the Guinier and the Porod regimes, where is verified 
/ oc q~ Df [HI]. Various works can be found on the litera- 
ture which use light scattering: spherical vesicles contain- 
ing magnetic particles (magnetic liposomes) which gives 
a value of £>/ = 1.1 ±0.1 [2g|, organic-solvent-based mag- 
netic fluids compressing disperse magnetite nanoparticles 
with Df — 1.22 [2711 . and magneto-rheological fluids with 
D f = 1.21±0.14[25| andD/ = 1.24±0.04 [2f|. By means 
of the expressions explained in last section, it may be pos- 
sible to compare perimeter-based fractal dimension with 
three-dimensional fractal dimensions obtained using light 
scattering (34[. Helgesen et al [2l[ reported variation of 
fractal dimension with the magnetic strength for a 2D 
system of 3.6 /im magnetic particles contained in a cell of 
5 /im thickness. The fractal dimension is calculated using 
video-microscopy and EcQ. In fact, their obtained frac- 
tal dimension must not be considered three-dimensional 
but two-dimensional, because of their experimental sys- 
tem is clearly two-dimensional. The magnetic interac- 
tion is measured using the factor K^d — m 2 /8a 3 ksT. 
Comparing this quantity with our equivalent A factor, 
we obtain Kdd = A/2 using S.I. units. In their study, 
magnetic particles have a permanent magnetic moment, 
thus comparison between our experiments and theirs is 
complicated. In one case, they apply a external mag- 
netic field of H = 1 Oe, generating chains aligned with 
the external field. This experiment, the more compatible 
with ours, a value of D2 — 1.05 ± 0.03 for Kdd = 1360 
(A = 680) is obtained. Simulations have been reported 
on magnetic fluids including van der Waals, electrostatic, 
magnetic dipole and hydrodynamic inter-particle forces 
f32| . In that work, Chin et al obtain values for fractal 
dimension inside the interval Df = 1.1 — 1.3 for magnetic 



fields 0.1, 1 and 2 T. However, tangential magnetic dipole 
forces have not been included, so aggregates do not have 
a chain structure, nor are aligned with the magnetic field 
and therefore comparison with real experimental systems 
may be not applicable. Changes and temporal evolution 
of fractal dimension during agg regation have been stud- 
ied on latex micro-spheres [14j . and lately on magneto- 
rheological fluids [26[ . 

In Fig. [5J we display an example of the values of D2 
and Dp obtained along the aggregation process with the 
requirement of r > 0.98. We also represent several aver- 
age values of Db(2D) for three different experiments. A 
temporal dependence in projected fractal dimensions is 
not detected and a relatively uniform value for fractal di- 
mensions which do not depend of aggregation dynamics 
is found, i.e., the fractal dimension is constant in time, al- 
though dynamical quantities such as number of clusters, 
mean cluster size and others similar magnitudes follow a 
power law behaviour in this kind of systems [4(| Hgj . We 
do not observe temporal variation of these fractal dimen- 
sions during initial stages of aggregation. This is because 
free particles are not detected suitably and this reduce 
considerably the correlation factor on linear regressions. 
As we justify later in this work, we can interpret capacity 
dimension as a measure of the degree of cluster space fill- 
ing, while the perimeter-based dimension is linked with 
the clusters border roughness. When clusters become 
linear chains, these kind of aggregates do not change the 
way they fill the space (they are almost linear objects) 
and the box counting dimension do not vary. In the same 
way, chains' border roughness does not suffer special al- 
terations during aggregation and consequently D p does 
not show any variation. 

We calculate the projected fractal dimensions varying 
A and §id ■ The obtained values for D2 and D p are sum- 
marized on Table Q] on columns 3 and 4, the associated 
A and 4>2D values can be consulted on columns 1 and 
2. No dependency with A, 4>2D or with the ratio Ri/Rq 
is observed on these fractal dimensions. Therefore, we 
assume that making average values using all the experi- 
ments is a correct procedure. We also calculate the one 
dimensional dimension D\ for the chains, also with no 
variation or dependencies during aggregation, obtaining 
an average value of (D{) — 1.01 ± 0.03. One-dimensional 
fractal dimension gives some information about the irreg- 
ularity on the perimeter of the aggregate. For Euclidean 
dimension, we must have D\ = 1. In our case, we obtain 
a very similar value to this Euclidean result, probably 
because of the linear and field-aligned structure of the 
chains. Perhaps, as proposed Lee et al [34|, could be 
argued the existence of a relation among circularity C 
and Di, since C also measures the aggregate morphol- 
ogy as D\ does. Nevertheless, no correlation is observed 
between D\ and C for different kind of clusters, as Lee 
et al obtained. In our case, where aggregates are linear 
chains, C varies in the aggregation process with a power 
law behavior during aggregation (as it can be seen on 
Fig. [2]), while D-y keeps constant. 
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FIG. 6: Dependency of D B (2D) with the ratio R 1 /R dur- 
ing aggregation (Table [I] columns 6 and 3) . The continuous 
line represent a linear regression of the experimental data giv- 
ing D B (2D) =(-0.06 ± 0.01)Ri/Ro + (1.39 ± 0.02) with a 
correlation factor r =-0.873. 



Two-dimensional fractal dimensions Di and D p are 
more useful that one-dimensional fractal dimension D\ 
in terms of characterizing the formed chains from our 
magneto-rheological system. These dimensions not only 
give information about the morphology of the aggregates, 
but they also supply mechanical strength information 
about the aggregate. For two-dimensional fractal dimen- 
sion £>2, an average value of (D2) = 1.09 ± 0.02 is ob- 
tained. This value is quite close to the expected D 2 = 1 
for a linear object in an Euclidean geometry and is com- 
patible with the previous value obtained by Helgesen et al 
of D2 = 1.05 ± 0.03. For the perimeter-based dimension 
D p , we obtain an average value of (D p ) — 1.84 ± 0.02. 
Both values, D2 and D p , must be in an interval between 
1 and 2. For D p = 1, we have a perfect spherical object, 
whereas for D p = 2, we obtain a linear object. In our 
case, we obtain a value close to 2, so the aggregates are 
linear objects. The main practical difference among D2 
and D p is that D p is obtained from magnitudes which are 
calculated more directly (area and perimeter) and that 
Z?2 (and Di) is calculated by means of Feret's diameter, 
i.e., Dp reflects better the roughness of the boundary of 
the aggregate. For this reason, a value for D p is obtained 
not so close to 2 than D2 to 1 . 

In contrast to D±, D2 and D p , the capacity dimension 
seems to show a relation with the external parameters, 
so an average value cannot be calculated. In fact, we find 
a linear dependence with the ratio Ri/Rq as can be seen 
on Fig. [SJ As we mentioned on the Introduction, the 
ratio Ri/Rq informs about whether the chains initially 
aggregate by magnetic interaction (i?i > Ro) or by dif- 
fusion (i?i < Ro). In all our experiments Ri/Rq > 1 is 
verified, therefore, from the beginning, the aggregation 
process is dominated by the magnetic interaction among 
particles. Moreover, this result shows that the shape of 
the chains depends on the intensity of this ratio. Visu- 
ally, we observe that the finer and longer formed chains 
at the stationary state correspond to higher values of the 




Tima (b) 

FIG. 7: Four experiments as examples of power-law behaviour 
of the average chain-border height fluctuation W(t) in the 
initial stages of aggregation (t < 100) until saturation. From 
top to bottom: A = 131 and </>2d = 0.071; A = 703 and 
<p2D = 0.106; A = 252 and fan = 0.145; A = 1688 and 
4>id = 0.088. 



ratio Ri/Rq, while wider chains are observed when the 
ratio is lower. Therefore, the shape of the clusters may 
depend of the intensity of this magnetic interaction be- 
tween particles. This implies that the capacity dimension 
measures how the chains fill the surrounding space. 

In the Introduction, we introduced various concepts 
related with the roughness to the clusters contour, in 
particular, the border height fluctuation W(t). For our 
experiments, we observe that W(t) is approximately con- 
stant during aggregation from 100 s, until the field is 
switched off (see Fig. [7]). In Table U column 7, we show 
time average values {W) t for each experiment in the satu- 
rated region. No dependence with A, <f>2D or Ri/Rq is ob- 
served. Therefore, as it happens with D p , the variations 
on the border of the clusters do not suffer high variations 
during aggregation and they seem not to depend on ex- 
ternal variables, with the exception of the initial stages 
of aggregation. The average value of this dispersion is 
0.21± 0.04 //m, approximately 1/5 of the particle diam- 
eter. This value can be considered high for linear chains, 
because a perfect spherical single particle should have a 
value of w close to 1/3 of the particle diameter. However, 
the quantities {W) t have been obtained making a clus- 
ter average in which the free particles and small chains 
count the same as the longest clusters. Similarly, the ini- 
tial value when no field is applied may vary as a function 
of the number of initial aggregates present in the system. 
Therefore, this average and experimental variations for 
ideal scenarios may reduce the value of W{t) for a per- 
fect chain composed of aligned spherical particles. We 
observe as the average chain-border height fluctuation 
W(t) follows a power-law behaviour on the initial stages 
of aggregation, as can be seen for various experiments on 
Fig. [JJ Consequently, we define the exponent /3 a , so dur- 
ing initial stages W(t) ~ t@ a is verified. We summarize 
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FIG. 8: Example of chains disaggregating for different times 
in the experiment A = 131 and 4>2D = 0.071. The contour 
around the chains is the border of the clusters detected by 
means of the image analysis, while the grey color is the origi- 
nal captured image. Left: t = s (when field is off), Middle: 
t = 20 s, and Right: t = 40 s 




P a obtained values on Table [T] on column 5. An aver- 
age value for this exponent is ((3 a ) = —0.03 ± 0.02. The 
high relative error of this magnitude is associated with 
the differences in number and morphology between the 
detected clusters in each experiment, specially regarding 
the free particles. 

In addition, and as a verification, we calculated the 
roughness factor a for chains with sufficient length. We 
obtain a average value of three aligned chains (~ 20 /im 
length) a = 0.11 in a case of A= 703 and <j>2D = 0.106 and 
for different times t = 4105.6 s, 5233.2 s and 5400 s. If 
we apply Ec.©, we obtain an average value of d = 1.88 
± 0.03, compatible with the value D p = 1.85 obtained 
for this experiment and with the average (D p ) = 1.84 
± 0.02. As we will see below, this method also gives, 
where it is applicable, concordance between d and D p in 
the disaggregation process. Therefore, we can interpreter 
perimeter-based fractal dimension as a measure of the 
roughness of the chains border. 

At this point we would like to compare our experimen- 
tal findings with previous results on magneto-rheological 
fluids using light scattering. However this comparison 
may not be possible because of the nature of the three- 
dimensional fractal dimension obtained by static light 
scattering and the fractal dimension obtained using pro- 
jected images as in our experiments. However, theoretical 
expressions showed in the Appendix has been developed 
to carried out such a comparison. An adequate knowl- 
edge on this subject would be interesting in colloids for 
understanding the dynamics of the aggregates and their 
morphological evolution. Therefore, with our results of 
£>2, we obtain an average value of (Df) = 1.419 ± 0.002 
using Ec. ljAlap . By means of Ec. (|Alb[ ), we obtain an 
average value (Df) = 1.58 ± 0.04. Both results seem 
to give overestimated values, taking into account that 
they are clearly linear aggregates, especially regarding 
Ec. (|Alb[ ). Using the linear regression result Ec. (|A2[) on 
the data from Sanchez et al, a more sounder value is 
obtained. An average gives (Df) = 1.28 ± 0.04. A re- 
sult quite close to this average is obtained by means of 
Ec. (|A3p . using a resolution of N p i X = 1072 px, a similar 
resolution to our experimental system. Therefore, we ob- 
tain an average three-dimensional capacity dimension of 



FIG. 9: D p during disaggregation when the field is switched 
off (t = 0) in three different experiments (circles: A = 131 
and (j>2D = 0.071; squares: A = 176 and <j>2D = 0.068 and 
triangles: A = 703 and <j>2D = 0.106). The little crosses cor- 
respond to the average chain-border height fluctuation W(t). 
The black filled points are d — 2 — a values for the experiment 
A = 703 and 4> 2D = 0.106 



(Db(3D)} — 1.34 ± 0.01. Then, the more similar value 
that we have obtained comparing with previous experi- 
ments using static light scattering, is (Df) = 1.28 ±0.04, 
by means of the linear regression fitted with the data 
by Sanchez et al Ec. (IA2|) . A little higher value is ob- 
tained with the average value for the capacity dimen- 
sion Db(3D), using Ec. (|A3[) , this may be because the 
box-counting method gives imprecise results, due to its 
dependency of the pixel resolution. 



B. Disaggregation. 

As it has been said before, it is difficult to obtain re- 
liable results on fractal dimensions, when we have free 
particles because of image resolution. That is why we 
cannot observe the variation on fractal dimensions while 
the chains are being formed by the aggregation of individ- 
ual particles. However, we can observe the disaggrega- 
tion process, i.e, the separation by Brownian motion and 
repulsive electrostatic forces of the particles that formed 
the chains when the applied field is switched off. Surpris- 
ingly, in the literature, this kind of processes have rarely 
been studied 0, |50| . Developed studies focus on ag- 
gregation, but not on the process of the aggregates that 
vanish and become free particles. This is very interest- 
ing for practical applications, for example, for designing 
an on-off dispositive based on magneto-rheological fluids. 
In Fig. [5J we show three snapshots of a group of chains 
during disaggregation at times t = (field is turned off) , 
t = 20 s and t = 40 s. The dark border along the clusters 
border is the contour extracted using our software. It 
can be seen how the free particles have appeared at 20 
seconds and how the clusters contour turns into itself. 

The evolution of D p in the disaggregation process for 
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three different experiments varying A and <p2D values is 
displayed on Fig. [5] The perimeter-based fractal dimen- 
sion D p changes from the value D p ~ 1.90 for linear 
chains to a minimum value D p ~ 1.65 after 20 seconds of 
Brownian motion. After that, a stabilization on D p (for 
t > 40 s) is observed. We have to say that this method 
presents several limitations at this point: first, the pres- 
ence of free particles reduces the correlation factor of the 
linear regression, making the result less reliable. Sec- 
ond, as it has been commented in the Introduction, this 
method could not be correct when the objects do not 
present anisotropy, as occurs when the chain disaggre- 
gates into free particles. Therefore, D p values may not 
be correct when t > 60 s during disaggregation. 

In Figj9l we also represent (little crosses) the evolu- 
tion of the average chain-border height fluctuation W{t). 
Two separated states can be observed on this magni- 
tude: the first, an steep growth from t = until t ~ 10 
s; the second, a soft decrease of the curve until satura- 
tion. Both behaviours present, separately, power-law be- 
haviours with, of course, different exponents. We respec- 
tively named these exponents as (3di and (3d2 and their 
values for each experiment have been summarized on Ta- 
ble U columns 9 and 10. Average calculations of both 
quantities give the following results: (Pdi) — 0.15 ± 0.06 
and (f3d2) = —0.11 ± 0.03. We interpret these two re- 
gions as follows: the first region represents the region 
where the rugosity of the clusters grows because of the 
simple movement of the particles inside the cluster. In 
this stage, the clusters maintain their individuality not 
breaking yet into little pieces. This behaviour means uni- 
vocalicity in the contour points. The second region ap- 
pears when the clusters begin to break into little clusters 
and the free particles appear. At this point the growth on 
the average cluster rugosity decreases, because the aver- 
age is calculated with the new clusters and free particles 
that appear during the process. A crucial point for cal- 
culating the roughness exponent a in our clusters using 
the height- height correlation function Ec.([7]) is the uni- 
vocalicity of the contour. That means that, each x point 
must be associated to one and only one height. When the 
disaggregation process is advanced and the free particles 
appear, the detected contour becomes not univocal and 
tends to twist itself up. When this behaviour occurs, it 
is not possible to use this method. For this reason, we 
observe that we can only calculate a values for t < 4 s. In 
Fig. [51 using black filled points, we plot some calculations 
of d = 2 — a (Ec.©) for a long chain in an experiment 
with A = 703 and 4>2D = 0.106. The points of d fit ad- 
equately to the represented D p points and values seem 
to be very similar. For that reason, we conclude that 
the perimeter-based dimension, D pi is a measure related 
with the rugosity of the clusters contour. 

Finally, we calculate some average values for the box- 
counting dimension during disaggregation. These calcu- 
lations can be seen on Fig. [TUJ As we see on the study of 
aggregation, the box-counting dimension is different for 
each experiment depending on the ratio R±/ Rq, therefore 
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FIG. 10: Box-counting dimension during disaggregation for 
three different experiments. Each case begins the process with 
different Db(2D) values, but the fractal dimension tends to 
a common temporal evolution while the chains lose its linear 
structure. 

the experiments begin the disaggregation process with 
different values of Db(2D). However, when the field is 
switched off, all the values tend to a single curve. This 
might be the expected behaviour, if we suppose that the 
capacity dimension measures how the clusters occupy the 
space, since the clusters tend to dissolve into little clus- 
ters and particles, regardless of the characteristics of the 
field and concentration at beginning of the aggregation 
process. 

V. CONCLUSIONS. 

In this work, we study experimentally the cluster mor- 
phology in a magneto-rheological fluid (MRF) under a 
constant uniaxial magnetic field. We calculate 2D fractal 
dimensions and study the contour roughness of the clus- 
ter using image analysis. We study two processes: the 
aggregation, where the rod-like chains formed by the ap- 
plication of the external magnetic field, and the disaggre- 
gation which occurs when the magnetic field is switched 
off. As far as we know, the disaggregation process has 
not been studied in detail in the literature, in spite of its 
interest on potential practical applications. 

During aggregation, we do not observe any regular 
temporal variation on fractal dimensions, because the 
method used in this work cannot be applied to the 
initial stages of aggregation. Using the area-perimeter 
method, we obtain the following average values for: one- 
dimensional fractal dimension, (D\) = 1.01 ± 0.03; two- 
dimensional fractal dimension, (D2) = 1-09 ± 0.02 and 
perimeter-based fractal dimension, (D p ) = 1.84 ± 0.02. 
We compare these values with previous experimental 
works. We show that the fractal dimension obtained by 
Helgesen et al [Uj is indeed a two-dimensional fractal 
dimension. We obtain a good agreement between their 
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TABLE I: 2D fractal dimensions and /3 exponents. (Units for 
(W) t are /j,m.) 
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value with a external applied field (£>2 = 1-05 ± 0.03) 
and our result for two-dimensional fractal dimension. We 
obtain that the box-counting dimension or capacity di- 
mension Db(2D) does not vary with time, but its value is 
different depending on the ratio R\j Rq between two char- 
acteristic lengths that measures a relative magnetic field 
strength for particle. We also calculate various quantities 
associated with the roughness of the contour or border 
of the clusters. We find that the average chain-border 
height fluctuation W(t) follows a power-law exponent, 
in the initial stages of the aggregation process, with an 
average exponent {(3 a ) = —0.03 ± 0.02. 



During the process of chain disaggregation, i.e., when the 
applied field is switched off and the chains vanish because 
of the Brownian motion and electrostatic repulsive forces, 
we observe that D p decreases approximately a 13% fac- 
tor when the field is switched off. Also, we show that 
W(t) presents two different power-law behaviours with 
exponents (J3 d i) = 0.15 ± 0.06 and (f3 d2 ) = -0.11 ± 0.03. 
Moreover, we show that the perimeter-based dimension 
is related with the fractal dimension d obtained directly 
by means of the roughness exponent a in Ec®. This 
means that we may interpret the perimeter-based dimen- 
sion as a measure of the clusters border roughness. Fi- 
nally, we show results on the evolution of box-counting 
dimension during disaggregation. The differences among 
experiments on this fractal dimension vanish as the clus- 
ters dissolve. This behaviour fits with an interpretation 
of Db (2D) as a measure of how the clusters fill the space. 

In conclusion, for a better understanding of temporal 
evolution and dependencies of fractal dimension in MRF 
and for confirming the results hereby exposed, it would 
be necessary to carry out more investigations. Specifi- 
cally, simulations of MRF could be implemented using a 
3D system with a real set of interactions including elec- 
trostatic, van der Waals, Brownian motion and complete 



magnetic dipole forces. Such a model could be quite 
useful for an adequate study of fractal dimensions and 
for a better knowledge of the relation between three- 
dimensional fractal dimension and projected fractal di- 
mensions at an anisotropic colloidal system like the used 
in this work. Moreover, we will focus further investiga- 
tions on temporal evolution of aggregates' box-counting 
and correlation dimension and its dependence with mag- 
netic strength and concentration for the initial stages of 
aggregation. 
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APPENDIX A: CALCULATING 
THREE-DIMENSIONAL FRACTAL DIMENSION. 

How the values of these 3D and 2D-projected fractal 
dimensions relate to each other is not an straightforward 
matter. As we explained before, Df is usually obtained 
directly by means of light scattering, which is a very ef- 
fective experimental method in the study of colloidal sus- 
pensions. However, it would be very useful to obtain this 
three-dimensional fractal dimension using image analy- 
sis techniques. Some elaborated techniques have been 
developed regarding this problem, for instance, the two- 
slopes dimensional analysis using cumulative size distri- 
bution [5l[ or the relative optical density method [H3 |. 
The main question at this point is how to correlate eas- 
ily these 2D projected fractal dimensions with the mag- 
nitude which is actually useful physically: the three di- 
mensional fractal dimension Dr. A possible answer could 
be that D pro j = min{2,D/} [l(, but this result has been 
obtained for mathematical infinite self-similarity fractal 
objects and might not be applicable for real finite frac- 
tals, such as aggregates. A mathematical approach to 
this argument can be consulted on (3f| . Some theoretical 
studies have been developed in the last years to explore 
the possibility of obtaining Df using the projected fractal 
dimensions. In these works, the authors attempt to find a 
relation between the three-dimensional fractal dimension 
and the projected ones, using simulations of aggregation 
of particles in a three dimensional space. Afterwards, 
they calculate Df and the projected dimension and try 
to find mathematical relations between them. An initial 
work of the above-mentioned corresponds to Jullien et 
al [53[, but their results may be not applicable to real 
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fractal aggregates [35|| . As far as we know, the most im- 
portant approaches to this question appear in the works 
by Lee and Kramer (34[ , Maggi and Winterwerp [35| and 
Sanchez et al [36J. In the first above-mentioned works, 
Lee and Kramer find the following expressions between 
D f with D 2 and D p [Hi]: 

D f = 1.391 + O.Ole 2164 ^ 2 (Ala) 
D f = -1.628L> P + 4.6 (Alb) 

Although these two expressions are praiseworthy, the re- 
gression coefficient for Ec. (jAlbj) should be higher than 
the obtained r = 0.83, concluding that a reliable value of 
Df by means of the projected dimension D p has been 
demonstrated. Moreover, Ec. (jAlbj) must verify the 
boundary conditions for Df and D p , i.e., for a spheri- 
cal object, D v — 1 and Df = 3, and for a linear object, 
Dp = 2 and Df = 1. The first boundary condition is 
approximately verified, but not the second. On the other 
hand, Ec. (jAlaj) . which correlates Df with D 2 , has a 
higher correlation factor (r = 0.96) and seems to be more 
reliable. 

The boundary conditions for Df and D p suggest a pos- 
sible approximate linear expression in the form Df ~ 
—2D p + 5, which is approximately obtained using the 
data provided by Sanchez et al from their astrophysical- 
related work. In their study, Sanchez et al calculate a 
relation between D p and Df using simulations of ran- 
dom aggregates aiming to obtain the fractal dimension of 
interstellar molecular clouds by analysing the projected 
images of these objects. They gather that these relations 
mostly depend on the image resolution N p i X , especially 
at low values of D p (1 < D p < 1.6). However, as it can 
be seen on their work, there is a intermediate zone where 
there is no apparent dependence on the image resolution 
for values of D p contained in the interval 1.6 < D p < 1.9. 
For these data, a linear regression can be calculated, ob- 
taining the following expression: 

D f = (-2.09 ± 0.13)L> P + (5.15 ± 0.22); (A2) 
for 1.6 < Dp < 1.9 

The Ec. (jA2p verifies the two above-mentioned boundary 



conditions between Df and D p and gives an interval of 
application independent to the resolution. A mathemati- 
cal approach for obtaining the three-dimensional capacity 
dimension Db(3D), from two-dimensional projections of 
fractal aggre gat es is present in the work by Maggi and 
Winterwerp |35[. For statistical self-similar and non- 
homogeneous fractals, as fractal aggregates, it is usual 
to study the generalized dimension of qth order, d q [55| . 
The fractal dimension associated with q = gives do, 
being this dimension the capacity dimension Db- Maggi 
and Winterwerp apply the current theory of projections 
of fractals for obtaining the three-dimensional capacity 
dimension, Db(3D), that depends on perimeter-based 
dimension in the following form: 



^ (3D) ^rab ; ,olD ' <2 (A3) 

where a(N p i X ) and b(N p i X ) are functions of the resolution 
Npi X . In the case of infinite resolution I — > oo, a = 9/8 
and b = 7/8 is obtained. Therefore, it is possible to 
apply Ec. (jA3j) with the image resolution of our exper- 
iment. For example, for a resolution of N p i X = 1024 
px, a — 1.483 and b — 1.035 are obtained. An impor- 
tant result deduced of this study is that it is not possible 
to obtain three-dimensional information from the two- 
dimensional capacity dimension Db(2D). This result 
limits the utility of the box-counting fractal dimension 
for being compared with another similar quantities. In 
previous sections, we explained our results for the cal- 
culation of 2D projected fractal dimensions by means of 
image analysis. Afterwards, we calculated and compared 
the three-dimensional fractal dimensions obtained by the 
formed chains using the expressions explained in this sec- 
tion. The result that seems to be more simple and ap- 
plicable to our system is the obtained by Sanchez et al 
[36j ] which provide an approximate linear expression in 
the form D p ~ —2D p + 5 on the range of 1.6 < D p < 1.9 
and which is independent of image resolution. This result 
provides a value for three-dimensional fractal dimension 
of (Df) — 1.28 ± 0.04, compatible with li ght scattering 
studies of similar frameworks on MRF (25l. |26|. 
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